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Abstract 

Seesaw mechanism provides a natural explanation of light neutrino masses through suppression 
of heavy seesaw scale. In inverse seesaw models the seesaw scale can be much lower than that in 
the usual seesaw models. If terms inducing seesaw masses are further induced by loop corrections, 
the seesaw scale can be lowered to be in the range probed by experiments at the LHC without fine 
tuning. In this paper we construct models in which inverse seesaw neutrino masses are generated 
at two loop level. These models also naturally have dark matter candidates. Although the recent 
data from XenonlOO put stringent constraint on the models, they can be consistent with data on 
neutrino masses, mixing, dark matter relic density and direct detection. These models also have 
some interesting experimental signatures for collider and flavor physics. 
PACS numbers: 14.60.Pq, 14.60. St, 95.35. +d, 14.80.Ec, 14.80.Fd 
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I. INTRODUCTION 



Seesaw mechanism is one of the popular mechanisms [lH3| beyond the standard model 



Mi 



er than 
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(SM) which can provide some explanations why neutrino masses are so much smal 
their charged lepton partner mass scales mp. In Type I and III seesaw models 
requires the existence of heavy right-handed neutrinos of a Majorana mass scale M. The 
light neutrino mass is of order mo(mii/M). The usual scale of the heavy right-handed 
neutrino mass M is expected to be super heavy which can be as high as the grand unification 
scale. It would be good if the seesaw mechanism can be tested by high energy colliders. 
The LHC can test theoretical models beyond the SM at an energy scale as high as 8 TeV 
at present and will reach 14 TeV in the future. If indeed the heavy seesaw scale is of grand 
unification scale, it is impossible to test seesaw mechanism directly at accessible collider 
energies. Theoretically it is interesting to see if the seesaw scale can be lowered to TeV 
range allowing direct probe of ATLAS and CMS experiments at the LHC. There are indeed 
special solutions which allow lower heavy right handed neutrinos of order TeV with the 
price of fine tuning of the parameters 4| . Although this is theoretically allowed, it loses the 
original motivation of naturally explanation for the lightness of neutrinos through seesaw 
mechanism. Radiative seesaw neutrino mass generation can easy the problem and at the 
same time provide the much desired candidate for dark matter when additional symmetry 
exists to stablize the dark matter candidate jl 6]. The inverse seesaw mechanism^] can 
also lower the seesaw scale. In the inverse seesaw model, the heavy Majorana neutrinos are 
replaced by the heavy Dirac particles. The light neutrino masses are of order //(m^/M) 2 . 
Here M is the heavy Dirac particle scale and fi is a Majorana mass of the heavy Dirac 
particles which are supposed to be small even compared with mo- It is clear that the 
heavy Dirac mass scale M can be much lower than that for the Majorna mass in the usual 
seesaw models naturally. If the inverse seesaw is also achieved by radiative correction, the 
heavy scale can be even lower js il]. There are also other mechanisms to further lower the 



scale by naturally having a small fi parameter, such as that discussed in Ref. [12j through 
extra warped dimension. Here we will study radiative inverse seesaw models. In achieving 
radiative mass generation, sometimes it involves introduction of new symmetries to forbid 
terms which may induce tree level neutrino masses. If the symmetry introduced is unbroken, 
there may be a stable new particle in the theory. This new particle may play the role of the 
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dark matter needed to explain about 23% of the energy budget of our universe [13]]. In this 
paper we study several simple models using a leptonic heavy Dirac multiple to facilitate 
radiative inverse seesaw neutrino mass and also to have dark matter candidate. 



II. TREE INVERSE SEESAW 

The inverse seesaw neutrino mass matrix M v is the mass matrix resulted from the effective 
Lagrangian 

L m = -u L m D N R - N L MN R - -N&rNr - ^N L ^ L N C L + h.c. (1) 

where vl is the light active neutrino, Nl,r are heavy neutrinos. 
In the bases (V£, N R , iV£) T , M v is given by 



' m D » 



(2) 



\ M n L J 

With the hierarchy , L ~ , R « m D « M, the hght neutrino mass matr.x m„, defined 
by L mass = -(l/2)v L m„vl, to order (m D /M) 2 is given by|8|] 

m u = m D M~ x ii L {M~ x ) T m^. (3) 

There are different ways to achieve inverse seesaw mechanism depending on where Nl }R 
comes from. We briefly outline two simple possibilities which may realize inverse seesaw at 
tree level. 

One of the simplest ways is to introduce right-handed N R and left-handed Nz singlet 
heavy neutrinos with a discrete N R — > N R , Nl — > —Nl Z2 symmetry and all other SM 
particles do not transform under this symmetry. The mjj term is generated through Yukawa 
coupling L L Y D HN R . Here H = (h+, (v H + h + ilf '/a/2 is the SM Higgs doublet which 
transform under the SM electroweak gauge group SU(2)i x U(l)y as (2, 1/2). H = ia2H* . 
vh is the vacuum expectation value (vev) of H. Ll = {yL-,^L) T '■ (2,-1/2) is the SM 
lepton doublet. The [il,r terms are from bare Majorana mass terms N R /j, R N R and Nl/alN^. 
Because the Z 2 symmetry, the bare Dirac mass term NlMN r is not allowed. In order to 
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generate a non-zero M term, one can introduce a singlet scale S transforming under Z 2 as 
S — > —S with a non-zero vev v s /y/2. In this case, the Yukawa term N L Y S SN R is allowed 
which generates a M given by Y s v s /y/2. 

One can also introduce a leptonic doublet D^ R : (2, —1/2) along with a singlet S and a 
triplet A : (3, -1) (A^ with A n = A , A 12 = A 21 = A~/y/2 and A 22 = A -- ) to realize 
the inverse seesaw. One can introduce a global C/(1)d symmetry to distinguish Dl and Ll- 
Under this symmetry D LR — > exp[iar>]D LR , S — > exp[— iar>]S, A — > exp[2iao] A, and other 
fields do not transform. We have the following Lagrangian relevant to neutrino masses 

L D = -L L Y D D R S - D L MD R - ^D L Y L D C L A - ^D c R Y R D R ti + h.c. (4) 

If both S and A develop non-zero vev's, the inverse seesaw mechanism is realized. This 
model, however, will have a Goldstone boson due to breaking of the global U(1)d symmetry 
which may be problematic. To avoid the existence of a Goldstone boson in the theory, 
extension is needed. Also in the above two models, no candidates for dark matter. 

In following sections, we will extend the two models discussed in this section to 
radiatively generate inverse seesaw neutrino masses. We will discuss the possibility of 
having dark matter candidates in these models. 

III. RADIATIVE TWO LOOP INVERSE SEESAW 

To avoid the appearance of massless Goldstone boson in the theory, a possible approach 
is not to allow the global symmetry to break and therefore no Goldstone boson emerges. 
Applying this idea to the model involving Dl R , S and A are then not allowed to have vev's. 
This however also firbids the light neutrinos to have non-zero masses at tree level. We have 
to extend the model. To this end, we introduce another singlet a which transforms under 
the U(1)d as a — > exp[2iaD]<J. We refer to this model as the U(1)d model. The allowed 
renormalizable terms in the potential Vd are given by 

V D = -pL 2 H H^H + X H (H^H) 2 + f4^S + Xs(S^S) 2 + fjL 2 a ^a + X a (^a) 2 

+ /^AtA + A^(AtAAtA) Q + ]T X^ij^j + (fi Sa S 2 a + X AaH HAa^H + h.c), (5) 

ij 
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where the sum J2a is over ai l possible i and j, and i to be one of the H, S, a and A. The 



/l 3 

allowed terms are: 



\P HA (H^HA^A)p + \ Ha {H^Ha ] a) + \ HS (H ] HS ] S) + Aa ct (A^A(tV) + A^WS) . (6) 



In the above the indices a and /3 indicate different ways of forming singlet. They are given 
by 

(A t AA t A)i = A^-AyA^A*, , (A t AA t A) 2 = A*.A^A£Ai 
(A*ArfH)i = A].^, ,111.11,, , (rtAH^H) 2 = A^-if*^ (7) 

In the above \x\ are all larger than zero. The potential only allows H to have a non-zero 
vev Vjj. The theory has an unbroken U(1)d global symmetry after spontaneous symmetry 
breaking from SU(2)l x U{1)y to U(l) em . At the tree level, light neutrinos are massless. 
Giving the above terms in and Vd, it is not possible to defined conserved lepton number. 
This is because that among the L^D^S, Dl.rD c l r A and S 2 a, and HAa^H vertices, there is 
always one vertex where lepton number is violated. For example, assigning Ll to have lepton 
number +1 and Dl,r to have X, if one demands conservation of lepton number as LlDrS, 
Dl,rD c l r A and S 2 a vertices, at the vertex HAa^H would violate lepton number by 2 units. 
One can demand other vertices to converse lepton number, but no matter what one chooses, 
the combination of terms proportional to Y D Y L n Su X^aHYD always violate lepton number by 
2 units. However, the global U(1)d symmetry is respected. Because of lepton number is 
violated, at loop levels, Majorana neutrino masses may be generated. We find that non-zero 
Majorana neutrino masses can be generated at two loop level shown in Fig.l. This two loop 
contribution violates the lepton number by 2 units. This two loop mass generation is similar 
to the Babu-Zee model two loop neutrino mass generation [3] but with the light charged 
leptons in the loop replaced by new heavy particles. The last two terms in the potential are 
crucial for light neutrino mass generation. 

We will now describe how to calculate the two-loop induced light neutrino mass. After 
H develops vev, mixing between A and a will be generated via the term HAa^H. The 
corresponding mass matrix for (A , a) T can be expressed as: 

Ml Mf 2 
Ml M| 2 
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FIG. 1: Two loop diagram for neutrino mass generation. 



with 



M 2 U = & + -X^vjj , M 2 2 =fi 2 a + -X Ha v 2 H , Ml = M 2 
One can diagonalize the mass matrix via: 



-Xa<thv h 



cos a sin a 1 


ii 




— sin a cos a 1 







with 



sin 2a 



Ml + M| 2 ± V /(M 1 2 1 -M| 2 ) 2 + 4M 1 2 2 M| 1 



2M 2 2 



V /(M 1 2 1 -M 2 2 2 )2 + 4M 1 2 2 M 2 2 1 
The light neutrino mass m v generated via two-loop diagram is given by: 

d A p d 4 q 1 1 1 



m,, 



Y^M 2 Y L fXsa cos a sin a 



1 



(2tt) 4 (2tt) 4 p 2 - m 2 g 2 - m 2 p 2 — M 2 q 2 — M 2 
1 



(8) 



(9) 



x 



(p - g) 2 - m\ (p - g) 2 



m 



(10) 
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The last factor in the above can be written as (m^ — m 2 2 )/((p — g) 2 — m 2 J((p— g) 2 — m 2 2 ). 



Using, sin a cos a (m? 



02/ 



^12 = \Xa<thv 2 Hi and neglecting the mass splitting between 



rrifa and m^ 2 in the denominator, we obtain 
XAa H Y L Y 2 fi Sa v 2 H M 2 f d 4 p d*q 



m,. 



2(M 2 - m 2 ) 2 



(2tt) 4 (2tt) 4 
1 1 



1 



1 



1 



[(p ~ q) 2 ~ m ^ >1 } 2 V 2 ~ m | Q 2 ~ m | P 2 ~ m | Q 2 ~ M 2 p 2 — M 2 q 2 — m? s 
1 1 . 

+ p2 _ M 2 g 2 _ M 2> ■ 

Our results for the two loop integral agree with that obtained in Ref. 



(11) 
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Carrying out the loop integrals, we finally obtain 

m » = 2(470^(1 - ml/M^ ^(m^ms, m s ) - g(m 01 , M, m s ) 



-g{m^ , m s , M) + g(m H , M, M)] , (12) 



where 



g(m 1 ,m 2 ,m 3 ) 



f 1 ^ 2 

l dx[l + Sp(l - /j 2 ) - - ^log/x 2 

Jo 1- fx 2 



with fi 2 = ax ^} x ^ , Q = — ^i- Sp(z) is the Spence function or the dilogarithm 

function defined as: 

Sp(z) = - f ln(1 ~ t} rft (13) 
Jo ^ 

To compare with the inverse seesaw mass formula, we rewrite the above in a matrix form 
in the bases where M is diagonalized 

VHY^ k (X AaH fi S ^ l )Y^V H 



~ TH «« > ( 14 ) 

m kk 



where K k \ is defined as: 



Kkl = Skl 2{Ly (1 - m 2 /M 2 ) 2 ^( m ^ m s, m s ) - g(m^ , M kk , m s ) 

-g(m^,m s ,M fefc ) + gira^ M kk , M kk )\ (15) 



If one identifies, effectively, mo = Y d vh, M = diag(Mu) and /i£ = with = 

{X^aH^sa)Y^ Kij, the light neutrino mass matrix is effectively an inverse seesaw mass form. 
We therefore refer this as radiative inverse seesaw mechanism. This model is different than 
those radiative inverse seesaw models discussed in Ref. 
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lOj where additional neutral heavy 



spin-half particles are introduced to generate radiative neutrino masses. 



16|. As 



The above formula can easily fit current data on neutrino mixing and masses 
an example, let us consider a simple case with Y L diagonal and Y D = UdUpmns- F° r 
the normal hierarchy, choose Yl = diag(l, 1.05, 2.01) x 10~ 2 , yr> = 10~ 2 , XavH = 0.1, 
li s „ = lOOGeV, rrifa = 300GeV, m s = 150GeV, M u = 500GeV, we can get all the three 
neutrino mass 2.804 x 10~ 2 eV, 2.936 x 10~ 2 eV, 5.636 x 10~ 2 eV, respectively. These are 
consistent with data. For inverted hierarchy case, we just need to replace Yl with Yl = 
diag(1.297, 1.317, 0.100) x 10 -2 , with all the other parameters unchanged, the neutrino 
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masses will be 4.90 x 10~ 2 eV, 4.98 x 10~ 2 eV, 3.78 x 10~ 3 eV, respectively. Again, these 
numbers are consistent with data. 

Along the same idea, the case with singlet heavy neutrinos discussed earlier, can also be 
modified to have two loop realization of inverse seesaw mechanism. To this end we impose 
on the theory a global U(l)s symmetry. The new particles beyond SM are: N^r : (1,0), 
r] : (2,-1/2), A : (3,-1) and S : (1,0). Under the U(l)s these particles transform as: 
Nl,r —> exp[iazs], r\ — > exp[— ias]r), A — > exp[—2ias), S — >■ exp[— 2ia]S. The Lagrangian Ls 
for the bare mass term and Yukawa couplings, and the potential Vg relevant for two loop 
neutrino mass generation are given by 

L s = -N L MN R - L L Y D N RV - l -N R Y R N R S - ~N L Y L N C L S^ + h.c. 

V s = HA V V& j V + XashHA^SH + h.c. + .., (16) 

where indicate other allowed terms. 

The light neutrino mass matrix can be obtained by replacing /is CT by /iAr?> an d by 
^ash in eqJIU In this model terms proportional to YdYlHa^ashYd violates lepton number 
by 2 units, for the same reasons for the U(1)d model. 

As long as the radiative generation of inverse seesaw neutrino masses is concerned the 
above two models are very similar. However, when considering dark matter physics, these 
two models have different features. We proceed to discuss them in the following. 

IV. DARK MATTER CANDIDATE 

Since in both the U(1)d and U(l)s models, the global symmetries are not broken, there 
are stable particles which may play the role of dark matter. Which one of the new particles 
is the lightest one depends on the parameter space and therefore determines which one plays 
the role of dark matter. 

In the U(1)d model, the heavy fermion particles have non-zero hypercharge and cannot 
play the role of dark matter. This is because that although dark matter relic density can 
be produced by dark matter annihilate into gauge particle with known interaction strength 
with sufficiently large dark matter mass, the direct detection rate from t-channel Z boson 
exchange would be too large. This possibility is therefore ruled out. The neutral components 
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of the scalar fields in the models are other possibilities which may be identified as dark 
matter. The neutral component A , has problem to play the role of dark matter due to its 
non-zero hypercharge. If the real and imaginary parts of the A masses m r and have a 
splitting 5 = m r — m.j, the non-zero hypercharge problem can be resolved by invoking the 
inelastic dark matter mechanism 17| , namely the scattering of a dark matter off nucleon is 



kinematically forbidden if the mass splitting 5 is larger than 100 KeV or so. In the U(1)d 
model, however, we find that it is not possible to generate a non-zero 5 for the real and 
imaginary parts in A , the inelastic dark matter mechanism is ineffective. 

The natural dark matter field is S. It does not have a non-zero hypercharge and does not 
mix with any particle having hypercharge. As long as dark matter properties are concerned, 
this model is very similar to the real singlet (darkon) model 18J and therefore similar dark 



matter properties 



19] and identical to the complex scalar singlet model 20j with degenerate 



mass for the real and imaginary parts of S. This is a typical Higgs portal model. Dark 
matter annihilation and detection are all mediated by Higgs boson. 

The term important is S^SH^H. Removing the would-be Goldstone bosons in H, we 
have 

XshS^SH^H = ^\sh(v 2 h + 2v H h + hh)SS^ . (17) 

The first term will modify the mass of S from /x| to Mf, = //| + XsHvfj/2. As long as 
dark matter annihilation and detection are concerned, the free parameters are: Md, X$h 
and also the Higgs boson mass. In the model, the Higgs boson h properties, its mass and 
its couplings to SM particles (fermions and gauge bosons), are very close to the SM Higgs 
boson hsM- The recent LHC data indicate that the mass is about 125 GeV2l|]. We will 
analyze the model using Higgs mass of rrih = 125 GeV. The dark matter relic density and 
direct detection constraints on the coupling Xsh and Md are shown in Fig. 2. Since now we 
have two degenerate components as dark matter, the constraint on Xsh from relic density 
is 1/ a/2 times smaller than the darkon model 



search from Xenon 100 



18]. The recent data on direct dark matter 



23j put the most strigent constraint on the allowed range for dark 



matter mass. The range of a few tens of GeV for dark matter mass is in trouble. However, 
dark matter mass about half of the Higgs mass and larger than 130 GeV is still allowed. 

The a field is also a possibility for dark matter since it does not have a hypercharge 
neither. It mixes with A after H develops vev through the term: X^nHAa^H . The 



9 



1.(100 
0.500 



0.100 
5 0.050 



0.010 
0.005 









P 1 ' 



100 150 
M if (GeV) 





^ DAMA 














CDMS 




"- _ ^XENONIOQ \ 








]25GcV 




XENON'O- — 











10 15 20 30 50 70 100 150 

M D (GeV) 



FIG. 2: Constraints on the coupling and dark matter mass from dark matter relic density and 



direct detection 
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231 ] for S as the dark matter with Higgs mass set to be 125 GeV. The projected 



XenonlT sensitivity is also drawn. 



lighter of physical particle which may play the role of dark matter will also has a non-zero 
Z coupling. However, in this case there is the mixing parameter to tune to satisfy the 
constraint. We find that as long as the parameter sin a is less than 10 -3 , the large direct 
detection cross section can be solved. We also checked that sin a of order 10 -3 can be made 
compatible with the neutrino mass generation requirement. With a < 10 -3 , the dark matter 
is dominated by the component a. The dark matter properties are similar to S. 

We now briefly discuss dark matter properties in the U(l)s model. In this model, the 
neutral scalars in rj, A and S, and the N are possible candidates for dark matter. The neutral 
components in rj and A have hypercharges and also there are no mass splitting between the 
real and imaginary parts, they have too large cross section for direct dark matter detection 
after fitting relic density requirement. The S field although does not have a hypercharge, 
it mixes with A , some fine tuning is needed to be compatible with direct dark matter 
detection data. The situation is similar to the case of a as the dark matter in the U(l)u 
model. This is similar to the case of a as dark matter as in the U(1)d model. 

The N field does not have hypercharge and may also play the role of dark matter. In 
this case, the dark matter relic density is achieved by t-channel exchange of rj induced iV N 
pair annihilate into lepton pairs, / + l~ and Vl vl- The annihilation rate is governed by the 
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FIG. 3: Constraints on the coupling and dark matter mass from dark matter relic density and 
direct detection for TV as the dark matter. The different curves are for rj mass to be 200 GeV, 300 
GeV and 500 GeV, respectively. 



Yukawa coupling Yd, the mass of 77 and also the dark matter mass Mr, = M. We have 
checked that there are parameter space where the correct relic density can be produced. 
In Fig.3 we show some correlations of the parameters which can produce the correct relic 
density. At the tree level, N does not couple to quarks. However, at one loop level, with L L 
and 77 in the loop N-N-Z coupling can be generated which can lead to sizeable dark matter 
direct detection cross section and at the same time satisfy the dark matter relic density 
constraint. The results are shown in Fig. 3. We again see that the recent XenonlOO 23] 
data put stringent constraint on the allowed range for dark matter mass. But the N can 
still play the role of dark matter with appropriate masses. 



V. CONCLUSIONS 

We have proposed two models, the U{1)d and U{l)s models, in which neutrino masses 
are generated through inverse seesaw mechanism at two loop level. In these models, a global 
U(l) is unbroken leading to a stable beyond SM new particle in each model. These stable 
new particles are natural candidates for dark matter. We find that these models can satisfy 
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current experimental constraints from neutrino masses, mixing, dark matter relic density 
and direct detections. 

Because in these models the neutrino masses are generated at two loop level and also 
inverse seesaw type, the seesaw scale can be as low as a few hundred GeV. This can lead to 
observable signatures. Before closing, we would like to make a few comments about some 
phenomenological implications of the models. 

One of them is related to Higgs properties. Although the Higgs couplings to SM particles 
are not modified at tree level, there are noticeable corrections at one loop level in the above 
two models. An important example is the modification to h — >■ 77. Because the existence of 
the two terms (A^AH^H) a in both models and (rjjrjiHjHj ,T]Jr]jHjHi) terms in the U(l)s 
model, terms like A ++ A h and A + A~h (ri + i]~h) will be generated after H develops vev, 
the h — > 77 can be modified. At present the experimental value 2j| for this channel is 
1.9 ± 0.5 (ATLAS) (1.56 ± 0.43(CMS)) times that predicted by the SM. The central value 
is higher than the SM prediction. With large enough A^ A , one may bring the value to close 
to the data. 

Another is related to probing the new degrees of freedom in the models at the LHC. In 
both models there are new charged particles, the A, rj and D fields. The particles in this 
multiplet can be pair produced via electromagnetic and weak interactions. However, in both 
models there are unbroken U(l) symmetries, the new particles cannot decay into pure SM 
final state making detection difficult. A possible signature is that the charged new particle 
decays into an SM particle and a dark matter. The SM particle is detected, but the dark 
matter carries away large transverse missing momentum and energy. For example for the 
U(1)d model, with S been the dark matter, can decay into a charged lepton / ± and the 
dark matter S. In the U(l)s model with N being the dark matter, r/ ± can decay into a 
charged lepton and the dark matter. 

Finally, there are potentially large FCNC effects in leptonic sector in these models. This 

is because that the Yukawa couplings Yd in both models can be of order O(0. 1), at loop level 

exchange D and S, and, N and S in the U(1)d and U(l)s models, respectively, can generate 

flavor changing radiative decay of charged lepton I — >■ l ,r y with branching ratios close to the 

n 

current experimental bound[6|. Also possible large // — > e conversion. Near future improved 
experiments can test these models^]. Detailed analysis will be presented else where. 
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